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Abstract 

In many mathematical models for pattern formation, a regular hexagonal pattern 
is stable in an infinite region. However, laboratory and numerical experiments are 
carried out in finite domains, and this imposes certain constraints on the possible 
patterns. In finite rectangular domains, it is shown that a regular hexagonal pattern 
cannot occur if the aspect ratio is rational. 

In practice, it is found experimentally that in a rectangular region, patterns of 
irregular hexagons are often observed. This work analyses the geometry and dynam- 
ics of irregular hexagonal patterns. These patterns occur in two different symmetry 
types, either with a reflection symmetry, involving two wavenumbers, or without 
symmetry, involving three different wavenumbers. 

The relevant amplitude equations are studied to investigate the detailed bifur- 
cation structure in each case. It is shown that hexagonal patterns can bifurcate 
subcritically either from the trivial solution or from a pattern of rolls. 

Numerical simulations of a model partial differential equation are also presented 
to illustrate the behaviour. 
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1 Introduction 



Pattern formation is a topic of intensive current research; a recent volume of 
this journal has been dedicated to the subject [14]. There are many experi- 
mental examples of pattern formation, including Rayleigh-Benard convection 
[6], directional solidification [8] and the Faraday experiment [11]. It is now ap- 
preciated that the amplitude equations describing the evolution and dynamics 
of patterns are generic, so that many different physical systems are governed 
by the same equations. These amplitude equations and their possible solutions 
are determined more by the symmetries of the problem than the details of the 
physics [3]. 
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Theoretical work on pattern formation has generally been concerned with re- 
gions of infinite horizontal extent, with attention focussed on modes with a 
single wavenumber. Analysis of the competition between rolls and hexagons 
shows that in systems without an up-down symmetry, hexagons are generally 
preferred. For hexagons, a sign change in the pattern gives a genuinely differ- 
ent solution, while for rolls a sign change is equivalent to a spatial translation 
of half a wavelength. This means that hexagons can occur at a transcriti- 
cal bifurcation but rolls can only appear at a pitchfork bifurcation. In systems 
with an up-down symmetry, only pitchfork bifurcations occur, and either rolls, 
squares, or hexagons may be stable, depending on the coefficients in the am- 
plitude equations. 

Recently, more complicated patterns involving four sets of rolls with the same 
wavenumber have been considered [5]. These share the symmetry properties 
of rolls and squares, so a sign change is equivalent to a spatial shift. 

Analysis of patterns involving a single wavenumber is not always successful in 
predicting patterns. Resonant mode interactions can be important near certain 
points in parameter space. These resonances generate additional quadratic 
terms in the amplitude equations and can lead to more exotic patterns with 
complicated dynamics [15]. 

In most examples of experimental interest, there is some up-down asymmetry 
which leads to a preference for hexagons at the onset of instability. How- 
ever, the finite geometry of the experimental system or numerical simulation 
does not in general allow a regular hexagonal pattern. What is observed in 
practice is a solution which is close to regular hexagons, but in fact involves 
slightly different wavenumbers. An example is the numerical experiments on 
compressible magnetoconvection by Weiss et al. (1996) [18], who found pat- 
terns of non-regular hexagons in square domains. Laboratory experiments on 
surface-tension-driven convection in square container also yield a variety of 
puzzling patterns [9]. 

Such irregular or non-equilateral hexagons can also occur in a region of in- 
finite horizontal extent with anisotropy, since breaking rotational invariance 
also breaks the symmetry of the hexagons. An example of this is convection 
influenced by a weak shear flow [7,1]. Another application is in the numerical 
modelling of hexagonal patterns [16]; if the numerical method does not pre- 
serve the symmetry of the hexagons then a qualitatively incorrect bifurcation 
diagram can be obtained. 

Theoretical work on non-equilateral hexagons has been carried out by Mal- 
omed, Nepomnyashchy and Nuz [12], who studied the governing amplitude 
equations in two cases. If the pattern involves two equal wavenumbers, a hexag- 
onal pattern can be stable in a region of parameter space near the onset of 
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instability, and can branch either directly from the trivial solution or from a 
roll pattern. In the case of three different wavenumbers, hexagons can only be 
stable at finite amplitude. The same amplitude equations have also been stud- 
ied in the context of convection in the presence of a mean flow [7,1], and for an 
anisotropic solidification problem [8]. However, none of these works has pro- 
vided complete bifurcation diagrams showing how the well-known picture for 
competition between rolls and regular hexagons is modified by the anisotropy. 

This paper describes the geometry of irregular hexagonal patterns in section 
2, considering the constraints imposed by a finite domain and relating the 
wavenumbers involved in the pattern to the number of hexagons seen. Sec- 
tion 3 considers the nonlinear dynamics of these patterns, concentrating on 
those aspects of the problem not covered by previous work. To illustrate the 
possible types of behaviour and to provide a check on the analytical work, nu- 
merical simulations of the modified Swift-Hohenberg model [17] are described 
in section 4. 



2 Geometry of hexagonal patterns 

This section describes the geometry of hexagonal patterns that can be ob- 
tained in a square box with sides of length L with either periodic or Neumann 
boundary conditions. 

For periodic boundary conditions, any pattern w(x,y) can be written as a 
Fourier series, 

oo oo 

w = J2A mn exp27ii(mx + ny)/L (1) 

m=l n=l 

and for Neumann boundary conditions = at x = 0, L; ^ = at y = 0, 
L) 

oo oo 

w— A mn {cosn{mx + ny)/L + cosii(mx — ny)/L). (2) 

m=l n=l 

Note that the possible patterns for Neumann boundary conditions are a subset 
of those for periodic boundary conditions in a box twice as large. This fact is 
often referred to as a 'hidden symmetry' [2,4]. Henceforth periodic boundary 
conditions will be assumed, since this includes the other case. The integers m 
and n in (1) will be referred to as 'wave integers'. The corresponding wavenum- 
ber is 

k = 2nVm 2 + n 2 /L. (3) 
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A hexagonal pattern is obtained by taking three modes whose wave integers 
sum to zero, i.e. 

m 1 +m 2 + m 3 = 0, rii + n 2 + n 3 = 0. (4) 

This is necessary if 'strong' resonance (i.e. a quadratic term) is to occur in the 
amplitude equations. 

A pattern of regular hexagons with three equal wavenumbers cannot occur in 
a square box. This result is not immediately obvious (in general the hexagons 
could be aligned at any angle to the box) but can be demonstrated as follows. 
If the three wavenumbers are equal then 

2,2 2,2 2,2 /r-\ 

m 1 + n 1 — m 2 + n 2 = m 3 + n 3 , (5) 

which can be rewritten using (4) as 

m \ + n \ — m l + n l — ( m i + m 2) 2 + (ni + n 2 ) 2 . (6) 

Expanding the brackets, these equations simplify to 

m\ + n\ = m\ + n\ = —2{mim 2 + nin 2 ), (7) 

showing that the sums of the squares of the wave integers are even. However 
it can be stipulated that not all the wave integers are even, since any common 
factors can be removed. There are two possibilities, each of which lead to a 
contradiction. If mi, rii are odd and m 2 , n 2 are even, then m\ + n\ = 2 mod 4 
but ml + n\ = 0mod4. Alternatively if mi, rii, m 2 , n 2 are all odd, then 
m\ + n\ = 2 mod 4 but 2{m x m 2 + nin 2 ) = mod 4. 

This result, that regular hexagons cannot occur in a square box, generalizes 
to the case of a rectangular box provided that the aspect ratio of the box is 
rational. In this case, if the box is of size Li in the x-direction and L 2 in the 
^/-direction and Li/L 2 = p/q where p and q are integers, the constraint of 
equal wavenumbers requires 

ml/L\ + n\/Ll = mi/Lj + n\lL\ = mj/Lj + nj/Ll (8) 

Multiplying through by L x L 2 pq, this becomes 

(miq) 2 + (nip) 2 = (m 2 q) 2 + (n 2 p) 2 = (m 3 q) 2 + (n 3 p) 2 (9) 

which is equivalent to (5) with m^ replaced by m^q and rii replaced by m t p, so 
that the problem has been reduced to the previous case. 
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Fig. 1. Hexagonal patterns with wave integers (a) (2, 1), (—2, 1), (0, —2); (b) (3,0), 
(-2,-2), (-1,2). 

Regular hexagons can however be obtained in a rectangular box if the ratio 
of the length of one side to the other is a multiple of v^; this choice has been 
used in numerical simulations to investigate the relative stability of rolls and 
hexagons [13]. 

Two examples of irregular hexagonal patterns are illustrated in Fig. 1. In each 
case the amplitudes of the three modes have been chosen to be real and equal, 
so that the function plotted is 

3 

w = ^2 cos(2n (rriiX + nit/)/ L). (10) 
%=\ 

Fig. 1 (a) shows the case where the wave integer pairs are (2, 1), (—2, 1) and 
(0, —2), in which four hexagons occur in the box, while the six-hexagon case 
with wave integers (3,0), (—2, —2), (—1,2) is shown in Fig. 1 (b). The former 
can occur with periodic or Neumann boundary conditions, but the latter can 
only occur with periodic boundary conditions. The six-hexagon pattern shown 
in Fig. 1 (b) has been found in numerical simulations of magnetoconvection 
in a compressible fluid by Weiss et al. (1996) [18]. 

It is apparent from these examples that there are two different possible sym- 
metries for these hexagonal patterns. The four-hexagon example, with two 
equal wavenumbers, has a reflection symmetry (x — > —x). The six-hexagon 
case, with three different wavenumbers, does not have a reflection symmetry. 
It is well known that symmetry has a crucial effect on the nonlinear dynamics 
of a system (e.g. Crawford and Knobloch 1991 [3]), and it will be shown in 
section 3 that the bifurcation structure associated with the two patterns is 
different. 
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There is a simple relationship between the number of hexagons N in the 
periodic box and the wave integers. This is 

N — | m\n 2 — m 2 ni\ . (11) 

From (4) it follows that any pair of the wave integers can be used in this 
formula. This result is demonstrated as follows: the planform function (10) is 
maximized when x and y obey 

m\X + riiy — pL, m 2 x + n 2 y = qL, (12) 

where p and q are integers. From (4) it follows that m 3 x + n^y is then also an 
integer multiple of L, so that w takes its maximum value of 3, corresponding 
to the centre of a hexagon. These equations give the lattice of points at which 
centres of hexagons occur. Three points on this lattice are 

(a) p = 0, q = 0, x = 0, y = 0; 

(b) p — 1, q — 0, x — n 2 L/ (mi^ — m 2 ni), y = -m2L/(min 2 — m 2 n\) ; 

(c) p — 0, q — 1, x — —n\Lj {m\U2 — m 2 ni), y = m\Lj{m\ni — m 2 ni) . 
These three points form a triangle that connects the centres of three hexagons. 
The area of this triangle is 

|(n 2 , —m 2 ) x (— ni,mi)| L 2 /2(m 1 n 2 - m 2 n 1 ) 2 = L 2 /2\m 1 n 2 - m 2 ni|,(13) 

using the formula for the area of a triangle in terms of the cross product of 
two vectors. This area is half the area associated with each hexagon, since 
each triangle connects three hexagons but each hexagon is connected to six 
triangles. Dividing the total area of the square L 2 by the area of each hexagon 
gives the formula (11) for the total number of hexagons in the box. 

To describe these hexagonal patterns it useful to have a parameter measuring 
how close the pattern is to regular hexagons. One such parameter is H = 
(^max ~ ^mirJAmax where /c max and k min are the maximum and minimum of 
the three wavenumbers; if H is small, the pattern is close to regular hexagons. 
Table 1 lists all hexagonal patterns in which H < 0.45 and iV < 25, giving the 
wave integers, the parameter H, the number of hexagons N and the symmetry 
for each. 



3 Nonlinear dynamics of hexagonal patterns 

This section describes the behaviour of the nonlinear amplitude equations 
governing the dynamics of irregular hexagonal patterns. 
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Table 1 

Wave integer combinations that generate hexagonal patterns. The parameter H 
measures the departure from regular hexagons and N is the number of hexagons 
in the periodic box. The final column indicates whether or not the pattern has a 
reflection symmetry. 



Consider a pattern-forming partial differential equation with a control param- 
eter r and a zero solution that exists for all r. Suppose that there is a marginal 
stability curve r c (k), on which the growth rate of perturbations to the zero 
solution vanishes, with a single minimum. There are many examples of such 
systems, including convection and the Swift-Hohenberg model [17]. Let r be 
the minimum value of r c (k), and let ko be the corresponding wavenumber. Now 
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for a finite box size L, only discrete wavenumbers given by (3) can occur, and 
in general these will not include k . For large values of L however, wavenum- 
bers close to k will be included, and the separation between wavenumbers in 
the vicinity of k scales as 1/L 2 . This allows a small parameter e to be intro- 
duced, representing the difference between fc and the three wavenumbers that 
generate the hexagonal pattern. With the scaling hi — ko + 0(e) for i — 1, 2, 3, 
the corresponding values of r c are r C j = r Q + 0(e 2 ), since tq is the minimum 
of r c {k). Writing r = r + e 2 r 2 , all growth rates are 0(e 2 ) so the appropriate 
time scale for the amplitude equations is T = e 2 t. 

It is now possible to write down the governing amplitude equations for the 
three Fourier modes obeying the resonance condition (4). Writing 

3 

w = ^2 Aj exp(27ii(rrijX + rijy)/L), (14) 
i 

it can be shown that the sum of the phases of the Aj tends to zero [12], so 
that with an appropriate choice of origin the amplitudes Aj can be taken to 
be real. The scaled amplitude equations are then 

A 1 = X 1 A 1 + A 2 A 3 , 

A 2 = \ 2 A 2 + A ± A 3 , (15) 

A 3 = X 3 A 3 +A 1 A 2 . 

Here, the dot represents the rate of change on the slow timescale T, the linear 
coefficients Ai, A 2 and A 3 are proportional to (r — r C j)/e 2 and the amplitudes 
Ai have been scaled by a factor e 2 . With these scalings, all terms in (15) are 
of the same order, and any cubic terms in the amplitude equations will be 
smaller by a factor e 2 . 

Note that the coefficients of the quadratic terms in (15) are equal. This is 
because (15) represents a small perturbation from the case of regular hexagons. 
The coefficients of the quadratic terms can be set to 1 by scaling the amplitudes 
Aj appropriately. 

The dynamics of the system (15) depends crucially on whether there are any 
additional symmetries. It is useful to begin by briefly reviewing the well-known 
case of regular hexagons, for which Ai = A 2 = A 3 . In this case the fixed points 
of (15) are the solution A 1 = A 2 = A 3 = 0, which is stable for Ai < 
and unstable for Ai > 0, and four equivalent hexagonal solutions A\ = ±Ai, 
A 2 = ±Ai, A 3 = —AiA 2 /Xi, for which the three eigenvalues are 2Ai, 2Ai and 
— Ai. The bifurcation at Ai = is transcritical, and the hexagonal solution is 
never stable. 
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In the case where Ai = A2, two of the wavenumbers are equal in magnitude 
and the resultant pattern has a reflection symmetry (e.g. Fig. la). The fixed 
points are 

(i) A 1 = A 2 = A 3 = 0. This is stable if Ai < and A 3 < 0, and undergoes 
stationary bifurcations at Ai = and A3 = 0. 

(ii) A 1 = ±a/AiA 3 , A 2 = ±VAiA3, A 3 = —AiA 2 /\ 3 . These four equivalent 
solutions only exist when A1A3 > 0. Its eigenvalues s obey s = 2Ai or s 2 — 
\ 3 s — 2A1A3 = 0. The product of roots of this quadratic is negative so this 
solution can never be stable. 

For the case where all three linear terms in (15) are different, the three 
wavenumbers are different and the hexagonal pattern does not have mirror 
symmetry (e.g. Fig. lb). The fixed points are 

(i) Ai = A 2 = A 3 = 0. This is stable if Ai < 0, A 2 < and A 3 < 0, with 
stationary bifurcations at Ai = 0, A2 = and A3 = 0, 

(ii) A\ = ±y/\ 2 \ 3 , A 2 = ±VAiA3, A 3 = — AiA 2 /\ 3 . Again there are four 
solutions of this type. This solution only exists when either Ai < 0, A2 < 0, 
A3 < or Ai > 0, A2 > 0, A3 > 0. The eigenvalues s obey s 3 — (Ai + A 2 + 
A 3 )s 2 + 4AiA 2 A 3 = 0. Since there is no linear term, the three eigenvalues obey 
s i s 2 + S 2 S 3 + s 3 s i = which shows that the solution can never be stable. 



3.1 Amplitude equations including cubic terms 

The analysis of the previous section, although asymptotically correct, has 
two drawbacks. Firstly, no stable nonlinear solutions are found, and secondly, 
solutions in the form of rolls (involving a single wavenumber only) are not 
found. These problems can be overcome by the addition of cubic terms to the 
amplitude equations. In general, this is inconsistent, since quadratic and cubic 
terms can only appear together if the amplitudes are 0(1), in which case terms 
of all order should appear in the equations. However, a consistent scaling can 
be obtained if the asymmetry in the problem (and hence the quadratic term) 
is small; this is a commonly used assumption [1,7,10]. The appropriate scaling 
is that the coefficient of the quadratic term is 0(e) (recall that e is the scale 
of the difference between the wavenumbers). The growth rates are 0(e 2 ) as 
before, and a consistent balance between linear, quadratic and cubic terms is 
obtained if the amplitude scaling is = 0(e). The scaled amplitude equations 
then become 

A 1 = AiAx + A 2 A 3 - A ± (Al + (3A 2 2 + (3A 2 3 ) 

A 2 = X 2 A 2 + A 1 A 3 -A 2 (A 2 2 + (3A 2 3 + (3A 2 1 ) (16) 
A 3 = X 3 A 3 + A X A 2 - A 3 (A\ + f3Al + f3A 2 2 ). 
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It can be assumed that the quadratic and cubic terms are equal in each equa- 
tion because the system (16) represents a small perturbation from the equa- 
tions for regular hexagons; any deviation appears at higher order. By choosing 
an appropriate scaling for both time and the amplitudes, the coefficients of 
the quadratic terms and the A? terms have been set to unity. The constant (3 
is problem-dependent and cannot be scaled out. For simplicity it is assumed 
that (3 > 1. This means that rolls are stable in the absence of the quadratic 
term. This is indeed the case for both Rayleigh-Benard convection and for the 
Swift-Hohenberg equation. 

It is helpful to consider first the familiar case of regular hexagons for which 
\i = X 2 = A 3 [6]. The bifurcation diagram for this case is shown in Fig. 2. 
Solutions in the form of rolls (e.g. A\ = ±\/Ai, A 2 = A 3 = 0) are stable for 
Ai > l/(/3— l) 2 . Hexagons (with A 1 = A 2 = A 3 ^ 0) appear at a transcritical 
bifurcation from the trivial solution where they are unstable, but gain stability 
though a saddle-node bifurcation at Ai = — 1/(4 + 8/3) = Asn and are stable 
for Asn < Ai < Atr = (/3 + 2)/(/3 — l) 2 . There are two regions of overlapping 
stable solutions and hence hysteresis; for Ai small and negative both the trivial 
solution and hexagons are stable, while for larger Ai both hexagons and rolls 
are stable. The rolls and hexagons are linked via a branch of mixed modes (e.g. 
Ai = A 2 7^ A 3 ) which are always unstable. The upper stability boundary of 
hexagons at Ai = Atr occurs at a transcritical bifurcation with D 3 symmetry, 
where the hexagons and mixed modes have a double zero eigenvalue; a centre 
manifold reduction near this point yields the normal form 

x = /ix + x 2 - y 2 , y = iiy-2xy (17) 

with the symmetries of rotation through 27r/3 and reflection y — > —y. This 
system is discussed in more detail later. 

The procedure for analysing (16) is similar to that for (15). The behaviour de- 
pends on whether or not any of Ai, A2 and A3 are equal, so these two cases are 
studied separately. Since the equations (16) have been considered by other au- 
thors [12,7], the details of the calculations are not given. Attention is focussed 
on points not described by the previous work. The case of symmetric patterns 
with two equal wavenumbers is described in section 3.2 and the asymmetric 
case with three different wavenumbers is discussed in section 3.3. 

3.2 Patterns with two equal wavenumbers 

The fixed points of (16) when Ai = A2 £1X6 clS follows. 

(i) A\ = A 2 = A3 = 0. This is stable if Ai < and A3 < 0, and undergoes 
stationary bifurcations at Ai = and A 3 = 0. 
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Fig. 2. Bifurcation diagram for regular hexagons. Solid lines indicate stable solutions 
and dashed lines represent unstable solutions. Pluses and minuses show the signs of 
the three eigenvalues of each branch. 

(ii) 'Ai-rolls' with A x = ±y/\^, A 2 = A 3 = 0; or A 2 = ±V>^, A 1 = A 3 = 0. 
These rolls are appear at a supercritical pitchfork bifurcation and are stable 
if A 1 >(l-A 3 + A 3 /3)//3(/3-l). 

(iii) 'A3-rolls' with Ai = A 2 = 0, A 3 = ^y/X^. These rolls also bifurcate 
supercritically, are stable if Ai < f3\ 3 — y/X^, and have bifurcations at Ai = 
f3\ 3 ± v^. 

(iv) A mixed mode which can be regarded as hexagons, for which A\ = ±A 2 . 
This solution also branches from the zero solution at a pitchfork bifurcation 
at Ai = 0, which is supercritical if 1 + (3 + I/A3 > 0. At onset, A 3 is much 
smaller than A\ and A 2 , so this pattern has a rectangular appearance. This 
solution can undergo saddle-node bifurcations, and branches from the A 3 -rolls 
at Ai = [3\ 3 ± y/X^. 

(v) A mixed mode in which all the amplitudes are different. This solution 
branches from the Ai-rolls at Ai = (1 — A3 + \ 3 (3)/(3(P — 1), and can never be 
stable. 

To draw the bifurcation diagram, the number of parameters can be reduced 
to one by supposing that (3 and the box size L are fixed but r is allowed to 
vary. This is equivalent to writing Ai = A3 + 5, with 5 fixed. There are two 
possible cases, according to the sign of 5, and these are shown in Fig. 3 and 
Fig. 4. In the bifurcation diagrams it is assumed that 5 is small, so that the 
bifurcations from the zero solution are close together. 

In the case 5 > (Fig. 3) the double bifurcation from the origin occurs first 
as Ai is increased. If 5 < 1/(1 + f3), a mixed mode with A 1 = ±A 2 bifur- 
cates subcritically from the origin and then gains stability in a saddle-node 
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Fig. 3. Bifurcation diagram for hexagons with mirror symmetry in the case where 
the double mode bifurcates first (5 > 0). The mixed modes are labelled according 
to whether all three amplitudes are different, or two are equal. 




Rolls 
Rolls 

Equal 

Equal 
Different 

Equal 



Fig. 4. Bifurcation diagram for hexagons with mirror symmetry in the case where 
the single mode bifurcates first (5 < 0). 

bifurcation. For 5 > 1/ (1 + /3), the mixed mode bifurcates supercritically. This 
mixed mode loses stability when it bifurcates to the (unstable) solution with 
all amplitudes different. Both roll solutions become stable at large Ai. 

For 5 < (Fig. 4), A3-rolls are stable at small amplitude but lose stability at 
a pitchfork bifurcation to the mixed mode with Ai = ±A 2 . This bifurcation is 
subcritical for small values of \5\. The mixed mode gains stability at a saddle- 
node bifurcation but then loses stability at a second saddle-node bifurcation 
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and rejoins the roll branch, at which point the rolls regain stability. For larger 
values of \S\ the region of stable mixed modes decreases and for 5 < — 1/4(/5 — 1) 
there are no stable mixed modes, so that rolls are always stable. 

It is of interest to note that these diagrams can be determined to a considerable 
extent just by considering the splitting of the hexagonal case (Fig. 2) induced 
by the symmetry-breaking. The single branch of rolls splits into two branches, 
and both hexagon branches become mixed modes with A-y = ±A 2 . The mixed 
mode with A\ = ±A 2 is unaltered, while that with Ay = ±A 3 becomes a 
mixed mode with all amplitudes different. It follows that in both the cases 
5 > and S < 0, at large Ai there must be two stable roll solutions, three 

mixed modes with A\ = ±A 2 (two with a +H — stability and one with H ), 

and one mixed mode with all amplitudes different. This approach can be used 
to construct Figs. 3 and 4 except near Ai = 0, where the amplitudes are small, 
and near Ai = A T r = (/? + 2)/(/3 — l) 2 . Near this latter point, where in the case 
of regular hexagons there is a transcritical bifurcation (17) with D% symmetry, 
the rotational symmetry is broken, but the reflection symmetry (y — > —y) is 
retained. A centre manifold reduction of (16) near Atr for small S leads to the 
normal form 

x = fix + x 2 - y 2 - 5, y = fiy-2xy, (18) 

where \i is proportional to Ai — Atr and x and y are rescaled forms of 2A 3 — 
A\ — A 2 and A\ — A 2 respectively. In the case of regular hexagons (5 = 0) 
three mixed modes cross through the hexagons as fi passes through zero and 
the phase portraits £1X6 clS shown in Fig. 5(a). For 5^0 the behaviour of 
(18) depends on the sign of 5. For S > 0, two pitchfork bifurcations occur at 

li = ±2-y/<5/3 and the sequence of phase portraits is as shown in Fig. 5(b). For 

5 < 0, two saddle-node bifurcations occur at /i — ±2v^— 6; the phase portraits 
are shown in Fig. 5(c). This analysis enables the correct connections to be 
made to the various branches in the Figs. (3) and (4). 

3.3 Patterns with three different wavenumbers 

In the case where the pattern involves three different wavenumbers, Ai, A2 and 
A 3 are all different. The fixed points of (16) are: 

(i) Ai = A 2 = A 3 = 0. This is stable if A x < 0, A 2 < and A 3 < 0, and 
undergoes stationary bifurcations at Ai = 0, A 2 = and A 3 = 0. 

(ii) Three branches of rolls, for example A\ = ±y/X[, A 2 = A 3 = 0. These are 
stable if Ai + A 2 < 2/3 Ai and (A 2 - /3Ai)(A 3 - f3\i) > Ai. 

(iii) A mixed mode in which the three amplitudes are all different, which 
bifurcates from the roll solutions. 
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Fig. 5. The splitting or 'unfolding' of the transcritical bifurcation with broken D3 
symmetry (18). Phase diagrams are shown with \x increasing from left to right, for 
the three cases: (a) 5 = 0: Three mixed modes cross through the hexagons, (b) 5 > 0: 
Two mixed modes with all amplitudes different disappear at a pitchfork bifurca- 
tion with a mixed mode with two amplitudes equal, then reappear at a pitchfork 
bifurcation with the other mixed mode, (c) 8 < 0: The two mixed modes with two 
amplitudes equal disappear and then reappear at two saddle-node bifurcations. 

With a fixed box size, we can set Ai = A2 + C* = A3 + 5 with a and 5 fixed. Since 
the three modes are equivalent, it can be assumed without loss of generality 
that 5 > a > 0, so that the rolls bifurcate in the order Ai, A 2 , A 3 as Ai 
increases. The following results can then be obtained for the roll solutions: 
The first branch of rolls to bifurcate is stable at onset. For small a and 5 
these rolls lose stability and then regain stability as Ai increases; for larger a 
and S these rolls are always stable. The second branch of rolls to bifurcate is 
unstable at onset and becomes stable at a bifurcation to a mixed mode. The 
third branch of rolls is also unstable at onset, has two bifurcations to mixed 
modes and is then stable. 

The bifurcation diagram for small a and 5 is shown in Fig. 6. There are up 
to five distinct branches of mixed modes. Three of these correspond to the 
mixed modes in Fig. 2 and the other two correspond to the hexagons. Near 
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Fig. 6. Bifurcation diagram for hexagons without mirror symmetry. 

the transcritical bifurcation at Ai = Atr a centre manifold reduction yields a 
second-order system analogous to (18): 



Here, both the rotation and reflection symmetries have been broken. The 
constants 7^ and 7^ are related to a and S but their value is not important. 
The behaviour of (19) is similar to that shown in Fig. 5(c) (but without the 
mirror symmetry), so that two of the mixed modes undergo a pair of saddle- 
node bifurcations, while two undergo no bifurcations, as shown in Fig. 6. 



4 Numerical simulations of the asymmetric Swift-Hohenberg equa- 
tion 

To illustrate the results of the preceding work, this section describes numerical 
simulations of the Swift-Hohenberg equation [17] modified by the addition of 
a quadratic term: 



This equation is probably the simplest with the required properties of a pre- 
ferred wavenumber (k = 1) and an asymmetry between u and —u, and (20) 
or variations of it are commonly used as models for pattern formation (e.g. 
[10]). A Fourier mode with wavenumber k has growth rate r — (1 — k 2 ) 2 . By 



x = fix + x 2 — y 2 — 7 



y = Hy-2xy- 7 W . 



(19) 




(20) 
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(a) 



(b) 




Fig. 7. Numerical solutions to (20) with s = 0.25. (a) r = 0.03, L = 23, showing 12 
irregular hexagons, (b) r = 0.05, L = 21.25, showing 9 irregular hexagons. 

substituting three Fourier modes of the form (14) into (20), the amplitude 
equations (16) are obtained with (5 — 2 and \ = 3(r — (1 — k 2 ) 2 )/As 2 . 

A Fourier spectral method was used to solve (20) with periodic boundary- 
conditions in the region < x, y < L. In Fourier space the linear parts of 
(20) can be solved exactly since the Fourier modes decouple. The nonlinear 
terms were integrated using the second-order Adams-Bashforth method. The 
initial condition chosen was a small-amplitude random perturbation to the 
zero solution. 

For the case of hexagons with reflection symmetry, the pattern with 12 hexagons 
with wave integers (4, 0), (—2, 3) and (—2, —3) was studied. In this case, there 
are two different possible bifurcation diagrams (Figs. 3 and 4) according to 
whether the double mode or the single mode bifurcates first. For L = 23, the 
(—2, 3) and (—2, —3) modes bifurcate first and therefore the appropriate bi- 
furcation diagram is Fig. 3. However for these parameters, 5 = 0.445, so from 
the analysis of section 3.2 we expect the bifurcation to irregular hexagons to 
be supercritical. The numerical results show 12 hexagons for 0.002 < r < 0.03 
and pure (3, 2) rolls for r > 0.05, with both solutions stable for r = 0.035 and 
r = 0.04. This is consistent with Fig. 3. The solution for r = 0.03 is shown in 
Fig. 7(a). 

For L = 24.5, the (4, 0) mode bifurcates first and the appropriate bifurcation 
diagram is Fig. 4. Rolls with wave integers (4, 0) were found for < r < 0.01 
and r > 0.025, and a hexagonal pattern including the (—2, 3) and (—2, —3) 
modes was found for 0.015 < r < 0.02. 
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Fig. 8. Number of hexagons N against box size L for the equation (20) with r = 0.02 
and s = 0.25. N = indicates that the zero solution is stable and N = 1 indicates 
that rolls are stable. The dashed line is the curve ./V = ^/3L 2 /8tt 2 , obtained by- 
dividing the area of the box by the area of a regular hexagon. 

The solution with 9 hexagons, involving the wave integers (3,2), (—3, 1) and 
(0, —3) was chosen to study the case without reflection symmetry. This so- 
lution was found when the box size L was 21.25, and at this value of L 
the (—3, 1) mode bifurcates first. A solution with 9 hexagons was found for 
0.02 < r < 0.055 and rolls with wave integers (3, 1) were found for r > 0.06. 
For r = 0.057 both solutions were found to be stable. Referring to the bifur- 
cation diagram (Fig. 4), rolls are also expected for small r; in fact however 
a hexagonal solution with 8 hexagons was found for r < 0.02. Interactions 
between hexagonal patterns with different numbers of hexagons are of course 
beyond the scope of the analysis of section 3. The solution for r = 0.05 is 
shown in Fig. 7(b). 

One question of interest which is difficult to approach analytically is which 
of the many possible hexagonal solutions listed in table 1 are observed as the 
box size L is increased. To investigate this question a sequence of runs was 
carried out with r = 0.02 and s = 0.25. For each value of L, the computation 
was continued until a stationary solution was obtained. By looking at the 
dominant Fourier modes the solution was then identified as either rolls or 
hexagons and the number of hexagons was determined using (11). Since there 
is the possibility of multiple stable states, the computation was repeated ten 
times for each value of L, with a different small-amplitude random initial 
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condition in each case. Values of L between 6 and 60 were used, with a step 
of 0.25 in L. The results are summarised in Fig. 8, which plots the number 
of hexagons against the box size. Where the number of hexagons is zero this 
indicates the zero solution, and a solution in the form of rolls is shown as one 
hexagon. Note that there are regions in which both rolls and hexagons are 
stable, and also regions in which two or even three different types of hexagons 
are stable. 

Of the hexagonal solutions in table 1, those with 4, 8, 9, 12, 15, 16, 20, 23 
hexagons were found. Note that, as might be expected, these are the solutions 
with low values of H, i.e. those that are closest to regular hexagons. There 
appears to be a preference for the symmetric hexagons. However this is no 
longer true at large values of the box size, where there are roughly as many 
asymmetric patterns as symmetric ones. 



5 Conclusions 

The main results of this paper are as follows. For numerical experiments 
with periodic boundary conditions, or laboratory experiments with Neumann 
boundary conditions, a perfectly regular hexagonal pattern is not permitted 
in a square box or any box with rational aspect ratio. Instead, patterns of 
irregular or 'non-equilateral' hexagons are observed. These patterns are com- 
posed primarily of three wavenumbers, and occur in two types. If two of the 
wavenumbers are equal then the resulting pattern has a mirror symmetry, but 
if all the wavenumbers are different it does not. There is a simple relation- 
ship (11) between the number of hexagonal cells in the pattern and the three 
wavenumbers. 

The nonlinear dynamics of these patterns is controlled by the amplitude equa- 
tions (16). The results of the analysis of these equations are summarised in 
the bifurcation diagrams of Figs. 3, 4 and 6. These provide a more complete 
picture of the dynamics than has been given by previous studies of the equa- 
tions [1,7,8,12]. A useful technique for clarifying the bifurcations that occur is 
the centre manifold reduction near the transcritical bifurcation with D 3 sym- 
metry. This reduces the algebra of the analysis considerably, enabling firm 
conclusions to be drawn regarding the connections between different solution 
branches. 

The analysis is closely related to that of Proctor and Matthews [15] who inves- 
tigated the interaction of modes with wave integers (0, 1), (1, 0) and (—1, —1), 
which give a pattern with the topology of hexagons (isolated upflows and con- 
nected downflows) but with square symmetry. Some of the results are similar, 
for example the fact that the quadratic resonances can lead to subcritical be- 
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haviour and a preference for three-dimensional patterns over two-dimensional 
rolls. 

This work was motivated by numerical simulations of compressible magne- 
toconvection [18], in which hexagonal patterns with six rising plumes were 
found. The results obtained suggest that for very small Rayleigh numbers, 
rolls should be stable, and that for larger Rayleigh numbers rolls should again 
become stable. However it must be borne in mind that the equations (16) 
are only valid if the up-down asymmetry (the departure from the Boussinesq 
approximation in the case of convection experiments) is small, and this is 
generally not the case. There is much more complicated dynamics in these 
numerical experiments and further analytical work is required to understand 
and interpret them fully. 

Acknowledgements 

I am grateful to the University of Nottingham, the Royal Society and the 
Nuffield Foundation for financial support. This work has benefited from dis- 
cussions with S. Cox, M. Proctor, A. Rucklidge and N. Weiss. 

References 

[I] S.M. Cox, IMA J. Appl. Mathematics (1997). 
[2] J.D. Crawford, Phys. Rev. Lett 67(4) 441-444. 

[3] J.D. Crawford and E. Knobloch, Ann. Rev. Fluid Mech 23 (1991) 341-387. 

[4] J.D. Crawford, J.P Gollub and D. Lane, Nonlinearity 6 (1993) 119-164. 

[5] B. Dionne, M.Silber and A.C. Skeldon, Nonlinearity, in press. 

[6] M. Golubitsky, J.W. Swift and E. Knobloch, Physica D 10 (1984) 249-276. 

[7] P. Hall and R.E. Kelly, Physical Review E 52 (1995) 3687-3696. 

[8] R.B. Hoyle, G.B. McFadden and S.H. Davis, Phil. Trans. A 354 (1996) 2915- 
2949. 

[9] E.L. Koschmieder and S.A. Prahl, J. Fluid Mech. 215 (1990) 571-583. 

[10] C. Kubstrup, H. Herrero and C. Perez-Garcia, Phys. Rev. E 54 (1996) 1560- 
1569. 

[II] A. Kudrolli and J.P. Gollub, Physica D 97 (1996) 133-154. 



19 



[12] B.A. Malomed, A. A. Nepomnyashchy and A.E. Nuz, Physica D 70 (1994) 357- 
369. 

[13] P.C. Matthews, M.R.E. Proctor and N.O. Weiss, J. Fluid Mech. 305 (1995) 
281-305. 

[14] Pattern Formation and Instabilities in Continuous Dissipative Systems, Physica 
D 97 (1996). 

[15] M.R.E. Proctor and P.C. Matthews, Physica D 97 (1996) 229-241. 

[16] A.C. Skeldon, K.A. Cliffe and D.S. Riley, J. Comp. Phys, in press. 

[17] J. Swift and P. Hohenberg, Phys. Rev. A15, (1977) 319-328. 

[18] N.O. Weiss, D.P Brownjohn, P.C. Matthews and M.R.E. Proctor, Monthly 
Notices of the Royal Astronomical Society 283 (1996) 1153-1164. 



20 



